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AN EXAMINATION OF WAKE OSCILLATOR MODELS 
FOR VORTEX-INDUCED VIBRATIONS 


INTRODUCTION 


Experimental results for vortex-induced vibrations of both pivoted and non-pivoted 
cylinders in a crossflow are presented in this report. Simplified models involving two ordinary 
differential equations are reviewed and shown to exhibit many of the features of the experimental 
results. Results of perturbation analysis of the simplified models are presented and used to 
explain some qualitative features, such as resonant frequency locking and hysteresis, in terms of 
bifiircations in the slow flow. 

In an experiment by Wei,^ a cylinder attached to a leaf spring is placed at the bottom of a 
flow tank (figure la). Since the leaf spring restricts the motion of the cylinder, the end attached 
to the spring remains fixed, and the cylinder can only pivot around that point at an angle 
transverse to the oncoming flowing water. Since the leaf spring also provides a restoring force, 
the cylinder tends to stay in the upright position if there are no external forces applied. As water 
flows past the cylinder, vortices are shed along the length of the cylinder on alternating sides, 
providing a periodic force on the cylinder. 

Figure 2 depicts data from the experiment by Wei. Note that for a range of flow speeds, 
the frequency of vortex shedding is the same as the cylinder oscillation frequency. This effect, 
called resonance lock-in, involves large-amplitude periodic motions of the cylinder near its 
natural frequency. Outside of this resonance region, the vortex-shedding frequency varies 
linearly with the flow velocity via the Strouhal relation: 



( 1 ) 


where S is the Strouhal number, V is the flow velocity, and D is the diameter of the cylinder. 

The features observed in Wei’s pivoted cylinder experiment are similar to features 
observed in experiments where the entire cylinder is allowed to move transversely to the flow; 
figures la and lb illustrate the pivoted and non-pivoted cylinder systems, respectively. Figure 3 
shows the experimental results from a study by Feng^ of a non-pivoted cylinder in a crossflow of 
air. 


The data in figures 2 and 3 share many common features. Although the details of the 
fluid flow (including, for example, the coherent vortices) are expected to be different for these 
two t 5 q)es of experiments, the phenomenon of resonance lock-in is common to both. One feature 
present in the data in figure 2 that is not present in the figure 3 data is a part of the lock-in region 
where the frequency varies linearly with flow speed (labeled II and S = 0.18). This is because 
Wei’s study uses water as the fluid, and Feng’s study uses air as the fluid. The difference 
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between these two experiments can be quantified by using the ratio of cylinder mass to the mass 
of the displaced fluid, which is an important quantity in the behavior of vortex-induced vibration 
experiments. This mass ratio is low for Wei’s study using water and high for Feng’s study using 
air. For more on the role of the mass ratio, see Govardhan and Williamson^ and Khalak and 
Williamson."^ For an overview of vortex-induced vibrations, see Blevins.^ 



Figure la. Setup for Pivoted Cylinder Experiments 

X 



Figure lb. Setup for Non-Pivoted Cylinder Experiments 
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Response Amplitude and Frequency 



U/FnD 


Figure 2. Experimental Data for a Pivoted Cylinder in Water (Data from Wef). 
(Note that both the amplitude of cylinder oscillation (lower portion ofplot) and 
frequency (upper portion of plot) are displayed as functions of flow speed.) 
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Figure 3. Experimental Data for a Non-Pivoted Cylinder in an Air Flow 
(taken from Harden and Currie^) 
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In the following section, a model will be analyzed that was originally derived to explain 
resonance lock-in in the non-pivoted experiment shown in figure lb. However, since pivoted 
(figure la) and non-pivoted (figure lb) experiments each exhibit resonance lock-in, the model 
may be used to explain some features of Wei’s data (figure 2). 


ANALYSIS OF WAKE OSCILLATOR MODELS 


Semi-empirical models for vortex-induced vibrations of cylinders have been widely 
studied. A review article by Parkinson® presents a good summary of these models. A van der 
Pol-type oscillator is commonly used to represent the time varying forces on the cylinder due to 
vortex shedding. The idea was first suggested by Bishop and Hassan.’ These types of models 
are known as wake oscillator models. Hartlen and Currie^ first proposed the following model in 
non-dimensionalized form: 


X + 2Cx + x = a co^y , (2) 

y + co^y-a(oy +—y^ = bx, (3) 

0 ) 

where x is the dimensionless cylinder displacement andy is a representative fluid property, for 
example, pressure or lift coefficient. The dots represent differentiation with respect to time t, and 
0 ) is proportional to the flow speed of the system, of which a, C y, a, and b are parameters. 
Here the cylinder and its elastic restraint are modeled by a damped linear oscillator (equation 

(2) ), and the periodic vortex shedding of the fluid is modeled by a limit cycle oscillator (equation 

(3) ). These two ordinary differential equations are assumed to be linearly coupled. Note that the 
fluid drag force on the rod i:|x| has been omitted. 

Figure 4 shows the results of numerically integrating equations (2) and (3) using the 
original parameter values (equation (4)) in the study by Hartlen and Currie. These parameter 
values correspond to the experimental setup in figure lb for the data in figure 3: 

C = 0.0015, a = 0.02, y=|, a = 0.002, and 6 = 0.4. (4) 

Figure 4 was obtained by numerically integrating equations (2) and (3) for a long time 
interval until most of the transients expired and the remaining steady-state motion was periodic. 
Comparing figure 4 with figures 2 and 3, it can be seen that some of the features observed in 
these experiments are reproduced by the Hartlen-Currie model. A large cylinder-oscillation 
amplitude resonance region occurs when the vortex-shedding frequency is near the natural 
frequency of the cylinder. Also, the frequency of oscillation in this region is nearly constant at a 
value close to the natural frequency of the cylinder. 
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Frec[uency (divided by 

natural frequency) X Amplitude 







Various methods have been used to analyze the dynamics of wake oscillator models. 
Poore and Al-Rawi^ have studied the Hartlen and Currie model in terms of Hopf bifurcations. 
Langford and Zhan*® performed a more detailed study on a similar wake oscillator model also 
using the Hopf bifurcation theory. Cor less and Parkinson ’ have studied the Hartlen and 
Currie model using multiple scales perturbation techniques around the lock-in region. Using a 
variety of methods, Berger has studied a similar model with an additional cubic damping 
term in equation (2) and a cubic x coupling term in equation (3). 

Figure 5 is an illustration of the amplitude and frequency response obtained from analysis 
of the Hartlen and Currie model using the original parameter values of equation (4). Features 
have been intentionally exaggerated to make them more visible. Arrows show jumps in 
amplitude and frequency occurring in response to sweeps of parameter to (proportional to flow 
velocity), resulting in hysteresis. These jumps occur because of changes in stability and 
associated bifurcations of periodic motions. Figure 5 was obtained by studying the Hartlen and 
Currie model using multiple scales analysis and bifurcation analysis of the resulting slow-flow 
equations; see the appendix for details. 




Figure 5. Amplitude and Frequency Response of Wake Oscillator Model (Horizontal axis co 
is proportional to flow velocity. Vertical arrows represent jumps in response to slow sweeps in 
flow velocity and result in hysteresis. Solid lines represent stable branches of periodic 
motions, and dotted lines represent unstable branches.) 
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Corless^^ has studied the possible response diagrams for different parameter values in the 
Hartlen and Currie model. This analysis uses the multiple scales method and is restricted to 
small parameter values near the 1:1 resonance. The analysis also is restricted to examination of 
periodic motions and does not examine quasi-periodic motions or the possibility of chaos in the 
system. 


CONCLUSIONS 


As shown in this report, simplified models of vortex-induced vibrations have been 
extensively examined in the literature. Nevertheless, there are two areas in which further 
research is needed: 

1. Although a great deal of attention has been paid to periodic motions, very little 
research has been focused on more complicated dynamical behaviors, such as quasi-periodic 
motions and chaos. 

2. Recent experiments have shown that an important parameter affecting the qualitative 
behavior of vortex-induced vibrations is the ratio of cylinder mass to displaced fluid mass. 
Although wake oscillator models have successfully described experimental results for large mass 
ratios, little has been done to model systems with small mass ratios. 
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APPENDIX 

RESULTS OF FREQUENCY AND AMPLITUDE RESPONSE 


The results of the amplitude and frequency response obtained from analysis of the 
Hartlen and Currie model using the Hopf bifurcation theory and perturbation methods are 
complementary. The Hopf bifurcation theory gives approximations that are valid close to the 
origin; perturbation methods give approximations that are valid for small values of 
a, C, y, a, b, and co~\ (near the 1:1 resonance of the system). The results from perturbation 
methods allow the amplitude and frequency response of a system to be determined for a given set 
of parameter values. Analysis has shown that the dynamics of the system involve Hopf 
biforcations from the origin where limit-cycle oscillations are created and destroyed. See Poore 
and Al-Rawi^ for more on the Hopf bifurcation analysis. 

In figure 5, points A and B correspond to Hopf bifurcations of the origin where limit 
cycles are created. Point C corresponds to a saddle-node bifurcation of cycles where the stability 
of the limit cycle changes along the branch. Point D corresponds to a torus bifurcation where a 
quasi-periodic motion is created. In the slow flow at point D, a subcritical Hopf bifurcation 
occurs where an imstable limit cycle (quasi-periodic motion in original equation) is created as co 
increases and the stability of the equilibria (limit cycle in original equation) changes. Note that 
between points C and D there are two possible stable solutions. The arrows show the jumps in a 
hysteretic loop for increasing and decreasing flow speed (represented by co ). For a more 
detailed bifurcation analysis, see Guckenheimer and Holmes.^^ 

A calculation similar to the one performed by Corless and Parkinson is presented below. 
To apply the method of multiple scales to the Hartlen and Currie model, a rescaling of equations 
(2) and (3) must first take place: 


n S a , 4C£ 

a = e, 2C = ep, « =—. b = Ygb, y- - j, 

Yq BFq 

and 

y = Ygy, x = ex 0)^—1 +eS, (A- 1 ) 

where e is considered a small perturbation parameter and ^ is a measure of detuning off the 1:1 
resonance. Substituting equations (A-1) into equations (2) and (3) and dropping the hats on 
terms gives 


and 


x + x = c(-;ffi: + qy)+ o{e^ ), 


y + y= e{-Sy + y l-^y^^ + bx)+o(e^). 


(A-2) 


(A-3) 
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Applying the method of multiple scales to equations (A-2) and (A-3) gives the following 
slow-flow equations: 


dR^ _ aRy sin(^!)) pR^ 

(A-4) 

dr] ~ 2 2 ’ 

dRy _ R'y Ry , 

(A-5) 

dr] 2 2 2 


d(p aRy cos (^) bR^ sin {(p) 6 
dX' 2Rx 2Ry 2 ’ 

(A-6) 


where x»R^ cos{^-6j, y»Ry cos{^-6y), (p = 6y-d^, ^ = (ot, and ri = et. 

For more on the method of multiple scales, see Nayfeh and Mook’’ and Rand and 
Armbruster.^* Equilibria in the slow flow correspond to periodic motions of the original system. 
Equations (A-4) - (A-6) can be combined to produce a single polynomial in Rx and the 
parameters a, b, p, and S. The right-hand sides of equations (A-4) and (A-5) are set to zero, 
and then solved for sin^ and cos^. These expressions for sin cp and cos (p are then substituted 
back into the right-hand side of equation (A-6) and sin^ (p + cos^ (p = \, producing two equations 
each involving Rx, Ry and the parameters a, b, p, and d. These two equations can then be 
combined, eliminating Ry to give the following polynomial: 

K^Rl + kX + kXx +Ko = 0 , (A-7) 

where 

K6 = by, 

K 4 = abf(4pd-2S-3ab), 

K 2 = cPiP^ - abd^ + Ip^S^ - +pd^- Sab p^S 

+ abpS + p^ - 2 p^+p^ + 3aX p). 


and 


ATo = -o' (PS^ - abpS -abd+p^-2p^+p + a^b\ 
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Equation (A-7) produces the curves of the amplitude plot in figure 5. Determining the 
stability of the curves involves linearizing equations (A-4) - (A-6) around the equilibria and 
examining the eigenvalues of the system. The saddle-node bifurcation at point C occurs when an 
equilibrium point has a zero eigenvalue. The Hopf bifurcation at point D occurs when an 
equilibrium point has a pair of pure imaginary eigenvalues. These calculations are 
straightforward but result in lengthy expressions that have been omitted. In a similar fashion, a 
polynomial for the lock-in frequency can also be calculated, producing the second plot in figure 
5. 
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